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ABSTRACT 



We investigate the detectability of Zeeman broadening in optical Stokes I spectra of slowly rotating sun-like stars. To this end, we 
apply the LTE spectral line inversion package SPINOR to very-high quality CES data and explore how fit quality depends on the 
average magnetic field, Bf. One-component (OC) and two-component (TC) models are adopted. In OC models, the entire surface is 
assumed to be magnetic. Under this assumption, we determine formal 3cr upper limits on the average magnetic field of 200 G for the 
Sun, and 150 G for 61 Vir (G6V). Evidence for an average magnetic field of ~ 500 G is found for 59 Vir (GOV), and of ~ 1000 G for 
HD 68456 (F6V). A distinction between magnetic and non-magnetic regions is made in TC models, while assuming a homogeneous 
distribution of both components. In our TC inversions of 59 Vir, we investigate three cases: both components have equal temperatures; 
warm magnetic regions; cool magnetic regions. Our TC model with equal temperatures does not yield significant improvement over 
OC inversions for 59 Vir. The resulting Bf values are consistent for both. Fit quality is significantly improved, however, by using 
two components of different temperatures. The inversions for 59 Vir that assume different temperatures for the two components yield 
results consistent with - 450 G at the formal 3cr confidence level. We thus find a model dependence of our analysis and demonstrate 
that the influence of an additional temperature component can dominate over the Zeeman broadening signature, at least in optical data. 
Previous comparable analyses that neglected effects due to multiple temperature components may be prone to the same ambiguities. 
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1. Introduction 

Cool stars display a range of phenomena that are typical of 
magnetic a ctivity, such as da rk spots present at the stellar sur- 
face (see Berdyugina 2005, for a review), enhanced chro- 
mospheric emission in, e.g., the C a II H and K lines (e.g ., 
IWilsonl 1 1994 iPizzolato et all 120031; lHall & Lockwoodl [2004), 



flares (e.g .. [F avata et al.l 2000; Maggio et al. 2000) and coronae 
(Schmitt 2001). Although these manifestations of stellar activity 
and hence of the underlying magnetic field are well detected for 
large samples of stars, measuring the magnetic field itself turns 
out to be difficult on most cool stars. The reason is the complex 
small-scale geometry of the field, that is amply evident on the 
Sun. The mixture of nearly equal amounts of magnetic flux of 
opposite polarities on the stellar disk leads to cancellation of the 
circular polarization signal due to the longitudinal Zeeman ef- 
fect, the most straightforward diagnostic of a stellar magnetic 
field. 

Rapid stellar rotation partly removes the degeneracy be- 
tween opposite magnetic polarities that lie sufficiently far apart 
on the stellar surface and thereby enables detection of a net 
(circular) polarization signal. Consequently, the magnetic field 
can be mapped on the stellar surface from the variation of the 
spectral shape and amplitude of this polarization signal over a 
stellar rotation, with the help of the Zeeman Doppler Imaging 



* Based on observations collected at the European Southern 
Observatory, La Silla 



(ZDI) technique (ISemefl 989). which has been successfully em- 
ployed to reconstruct the large-scale distribution ofjhe mag- 
netic flux on the surfaces of a number of stars dDonatil 12008; 
iDonati & Landstreetll2009l) . 

In spite of its success, ZDI does have some limitations. First, 
unless the field is homogeneously and unipolarly distributed on 
large scales, ZDI detects only some fraction of the true magnetic 
flux on the stellar surface. Second, the same is true for the mea- 
sured field strength and hence the magnetic energy density. In 
general, the field strength is underestimated and consequently 
the magnetic energy density is even more strongly underesti- 
mate d, since it is proport ional to the square of the field strength 
(see iReiners & Basri 2009). Third, ZDI is limited to compara- 
tively rapidly rotating stars. The slower a star rotates (or rather 
the smaller its v sin i) the lower the spatial resolution of the stel- 
lar surface (and of the magnetic field) that can be achieved. 

Therefore, the method proposed and applied to cool stars by 
Robinson ( 1986T); lRobinson et al.ldl980 ). of looking for changes 
in the line profile shape of the intensity (Stokes I) due to Zeeman 
splitting, remains an interesting technique to complement ZDI. 
In particular, it allows an estimate of the total magnetic flux as 
well as of the field strength and filling factor to be obtained (and 
hence also of the magnetic energy density). 

Zeem an broadening in stellar spectra (ZB) was first mea- 
sured bv IPrestonl (Il97l|) to determine t he av erage surface mag- 
netic fields of Ap stars. Robinson et al. ( 1980) then measured ZB 
in cool stars using his Fourier-based method which was taken 
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up and extended by other authors dGravll 19841 lMarcvl[T982h . It 
was soon realized, however, that the assumptions made in these 
Fourier techniques were too crude. Thus, a shift towards ever- 
improving forward radiative transfer calculations occurred. The 
deve l opment of these techniques was driven, am o ng others, by 
ISaail dl988lll990blll996l) : iBasri & Marcvl (1 19881): iLinskv et al.1 
(l!994l) : lRuedi et all(ll997h : lJohns-Krull & Valentil OOOOh . 

The displacement of a circularly polarized cr-component 
from line center can be expressed by: 



AA f 



4.67 x I0~ n g eS B A 2 



[A]. 



(1) 



For values typical of optical spectral lines, Eq. Q] demonstrates 
that even a comparatively strong field field along LOS is hard to 
detect. Bf = 1 kG would produce a Zeeman splitting of merely 
Ab - 17 mA for the cr-components of a normal Zeeman triplet of 
a hypothetical line at 6000 A with g e g = 1 .00, or Ab = 42 mA in 
the case of a more sensitive g e ff = 2.50 line. This should be 
compared to the line width of a typical stellar line profile broad- 
ened thermally and due to turbulence by 2-6kms~ 1 (i.e. 30- 
60 mA), as well as by rotation between and 5 km s _1 (i.e. up to 
100mA). In addition, the resolution element at even very high 
resolving power, R, is comparable in size, e.g. 5 A - 27 mA for 
fl = 220000. 

The present paper has a two-fold purpose. On the one hand, it 
is in the tradition of increasing the realism and accuracy of mag- 
netic field measurements on cool stars. It applies state of the art 
inversion techniques to determine the magnetic field. Also, from 
observations of the Sun it is clear that magnetic fields are as- 
sociated with significant brightenings or darkenings, which are 
generally neglected when measuring stellar magnetic fields. In 
the present paper we have carried out first simple computations 
looking at the effects that including such a diagnostic would 
have. We do not propose that this is an exhaustive or final study, 
but it can be considered as a first indicator. Finally, these tech- 
niques are also applied to spectral time series of an active star 
that has so far not been studied in this respect. 



2. Data 

2.1. Instruments and reduction 

Spectra of the stars investigated were obtained with the Coude 
Echelle Spectrometer (CES) at the 3.6m ESO telescope at La 
Silla, Chile. This fiber-fed spectrometer achieved a spectral re- 
solving power of approximately 220 000 using an image-slicer 
and projecting only part of an order onto the CCD, resulting in a 
small wavelength coverage of approx. 40 A. 

Spectra from two data sets were used in this work: data set 

A, with 5770 < A < 5810A, observed on 13 October 2000 
(solar spectrum) and 2 October 2001 (HD 68456), and data set 

B, 6137 < A < 6177 A, observed on 7 May 2005 (HD 68456, 
59 Vir). Sample spectra from both sets are shown in Fig. Q] We 
refer to iReiners & Schmittl d2003l hereafter R&S03) for more 
details, as data set A was originally taken for the project de- 
scribed in that publication. 

Where necessary, continuum normalization corrections were 
performed by division of linear slopes fitted to the original con- 
tinuum. Spectra were shifted to laboratory rest-wavelengths, b y 
cross-correlation using the FTS solar atlas (Kur ucz et alJ ll984). 



Line 


\ 


Ion 


lno (o* f*\ 


Xe 


geB 


Set 


1 Al 






TeVl 
LC v j 




A3: 


5778.46 


Fel 


-3.53 


2.59 


1.21 




5783.06 


CrI 


-0.57 (-0.50) 


3.32 


2.00 




5783.11 


CrI 


-1.32 (-1.82) 


3.56 


0.75 




JOUO. / J 


T 


n qq 


A ^1 


O Q9 
\J.yL 




5806.63 


Fe I 


-2. 19 (-2.09) 


4.91 


1.50 




5806.73 


Sell 


-3.02 


1.36 


1.67 


B3: 


6151.62 


Fe I 


-3.30 


2.18 


1.83 




6165.36 


Fel 


1 ^ i 


A 1 A 


l.UU 




6165.12 


Fel 


-2.76 


5.09 


0.67 




6173.34 


Fel 


-2.88 


2.22 


2.5 




6173.03 


Eu II 


-0.86 


1.32 


1.70 


B4: 


Same as B3 plus 










6162.18 


Cal 


-0.17 


1.90 


1.25 



Table 1: Spectral lines and their parameters used in inversions. 
Horizontal lines group line sets, double horizontal lines fur- 
ther distinguish the two data sets. Spectral lines are separated 
by small spacings indicating which blends belong together. 
Dominant line components are listed first, blends follow by as- 
cending wavelength. Brackets indicate the closest log (g*f*) (os- 
cillator strength) literature value, if values outside the literature 
range were adopted. 



2.2. Objects, spectral lines and atomic line information 

In this work, we present results for 4 objects ranging in spec- 
tral type from F6V to G6V, including sunlight reflected from 
Ganymede (previously and hereafter referred to as the solar 
spectrum). These stars were especially suitable for our analysis 
for multiple reasons. First, identification of spectral line blends 
was facilitated in the effective temperature range cov ered by 
these stars and LTE atmospheres were readily available dKuruczl 
1992). Second, the analysis of stars with relatively low projected 
rotational velocities, v sin i < 10 km s , allowed us to fully ex- 
ploit the high resolving power of the data and minimizes effects 
due to differential rotation. Third, measured stellar parameters 
such as effective temperature, T e g, surface gravity, log g and el- 
emental abundances, [Fe/H] in particular, were available as con- 
straints and could be directly input into the inversion procedure. 
Last, but not least, measured X-ray luminosities, Lx, and previ- 
ously reported magnetic field detections could be co nsidered as 
indicators for the presence of magnetic fields <HIIr|[l990b) and 
thus helped to constrain the sample. 

Within the small wavelength coverage of the CES spectra, 
only few relatively blend-free lines could be identified. Among 
these, we searched for combinations of lines that provide strong 
constraints on atmospheric properties. For instance, lines of dif- 
ferent excitation potentials, Xe, have different temperature sen- 
sitivities and often probe different heights (depending also on 
their equivalent widths). Therefore they constrain temperature 
stratification. In order to exploit the differential effect of Zeeman 
broadening, we aimed to cover the largest possible range in g e g 
as well as to use the line with the greatest possible g e g . 

Most spectral lines could be identified in the postscript ver- 
sion of the FTS solar atlas available online. The rema ining lines 
and b lends were identified by searching the V ienna ( Kunk a et al.l 
120001 hereafter VALD) and iKurucz & Belli d!995h atomic line 
databases. These databases further provided the atomic line data, 
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Fig. 1: Sample CES spectra illustrating data quality. Continua are shifted by ±0.3 from unity for visibility. Blue dash-dotted vertical 
lines indicate spectral lines used for analysis. Upper panel: data set A. Upper black spectrum: sunlight reflected from Jupiter's moon 
Ganymede; lower red spectrum: HD 68456 (F6V, v sin i « 10km s _1 ). Lower panel: data set B. Upper black spectrum: inactive and 
slow rotator 61 Vir (G6V); lower red spectrum: active and faster rotator 59 Vir (GOV, vsim » 6.5kms~ 1 ). 



such as level configurations and log (g*f*) values required for 
the forward radiative transfer calculations, cf. Sect.[3] Effective 
Lande values, g e g, were cal c ulated from the atomic level config- 
urations following Beckers (1969). T he calculated effective and 
empirical Lande factors determined bv lSolanki & Stenfld (ll985) 
were found to agree for all but one line: Fel at 6165.36 A, for 
which g emp = 0.69 and g e g = 1 .00. However, this line contains a 
blend at 6165.36 A, for which g eff = 0.67, cf. TableQ] 

For some lines, ranges of log (£*/*) values were found in 
VALD. We aimed to determine the best log (g*f) value for our 
selected lines by means of a by-hand gauging procedure using 
the solar spectrum for data set A and that of 6 1 Vir for data set 
B, while taking care to stay as close to the literature values as 
possible. This approach was adopted in order to avoid multiple 
additional free fit parameters, since one log {g*f*) for each blend 
would have been required. An automated approach would have 
been dominated by degeneracy between oscillator strengths and 
stellar model parameters. 

Table[T]lists all lines used and parameters adopted in our in- 
versions. Horizontal spacings help distinguish between blends 
contributing to the shape of a given spectral line and different 
lines. Col. 1 lists abbreviations used in this paper for the respec- 
tive line sets. Col. 2 contains wavelength in A, Col. 3 the corre- 
sponding ion in spectroscopic notation, i.e. Fe I is neutral iron. 
Col. 4 lists log(g*/*) adopted after gauging with the closest lit- 
erature value in brackets, if a value outside the VALD range was 
adopted. Cols. 5 and 6 contain excitation potential of the lower 
level of the transition, Xe, and effective Lande factor, g e ff, respec- 
tively. The range in Lande factor covered by line set A3 is thus 
Ageff = 1 .08, and Ag e s - 1 .50 for line sets B3 and B4, not count- 
ing blends. More importantly, the maximum g e ff is larger by 0.5 
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7.50 
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68456 
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2.90 


5.37 


7.19 


0.35 


59 Vir 


115383 


6.57 






7.78 


0.70 
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115617 
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7.55 
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Table 2: Logarithmic elemental abundances adopted on scale 
where hydrogen equals 12. 



in line sets B3 and B4 than in line set A3. Therefore, a stronger 
sensitivity to ZB in the latter two data sets can be expected. 

Collisional damping parameters are der ived internally by 
SPINOR from cross-section s publ is hed by Ans tee & O'Maral 
dl995l) : iBarklem & O'Maral d 1997b : iBarklem et al I d 19981) and 
are not listed in the table. 

Table|2] lists the elemental abundances a dopted. Solar abun- 
dances used were those determined by iGrevesse & Sauval 
(1998) and served as the absolute scale for relative stellar lit- 
erature values found in Valenti & Fischer ( 200 51 from hereon 
V&F05) (59 Vir & 61 Vir) and ICavrel de Strobel etail (1 19971) 
(HD 68456). 

Lines that correspond to elements for which no abundance 
measurements were found were usually not included in the in- 
versions. For the Cr I line in A3 as well as some blends, e.g. Eu II 
in B3, abundances were estimated from an iterative fitting proce- 
dure that took into account published average metallicities. This 
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procedure is potentially flawed, since general abundances are not 
necessarily correlated with individual abundances. However, this 
flaw does not challenge the work presented here, since Fe I lines 
were the most important ones for our analysis and literature val- 
ues for iron were available for all stars. The Cr I line was of great 
importance only for the solar case, for which measurements of 
the chromium abundance (and also that of scandium) do exist. 
The blend due to Eu II was only of secondary significance, given 
its small impact on line shape. It was not discarded, however, 
since it provided continuum points for (at least the blue wing of) 
the Zeeman-sensitive 6173.3 A line, for which the red wing was 
cut off due to the presence of unidentified blends. 

One well known Zeeman sensitive line was omitted from 
the analysis, Fell at 6 149.24 A. This line is very sensitive to 
magnetic fields due to its particular (doublet) Zeeman pattern, 
and displays c l ear Zeeman s plitting in spectra of Ap stars (e.g. 
Mathys 1990; jMathvs et al.lTl997h and solar penumbrae (e.g. 
Lites et al.ll99ll) . Unfortunately, this line is affected by blends in 
both line wings, as can be seen in Fig.[T]and we were not able to 
adequately reproduce these blended wings. While the blue wing 
blend is quite obvious, the red wing blend requires closer inspec- 
tion. This blend is (at least in part) due to Ti I at 6149.73 A and 
is hard to see due to rotational broadening, blending in with the 
smaller unidentified signatures seen in the slower rotator 61 Vir. 
This leads to the wing of the line just barely missing the contin- 
uum. Due to these two blended wings, the 6149.24AFeII line 
profile used for inversion had either a badly defined continuum 
(the deviation from the continuum is around 1%), or lacked both 
line wings (if they were both cut off as was done for the red wing 
of the 6173 A Fe II line) which meant not having any continuum 
points at all. Using either approach, we were not able to reach 
the same level of fit precision as presented for the other lines 
used in this analysis. Therefore, we chose to omit this line from 
our analysis. 

Our results depend on the assumed Signal-to-Noise ratio 
(S/N), since (S/N)" 1 = cr; and 

z — ' err 

where cr, is significance, i line index, j data point index for 
a given line, /y measured and I'., calculated intensity normal- 
ized to the continuum value. We determined cr, from the mean 
continuum standard deviation of the spectra. The typical S/N 
ratio thus determined was 400 - 700, corresponding to cr be- 
tween 0.0025 - 0.00143. Less statistical weight was assigned to 
the pressure-sensitive Cal line at 6162.18 Aby lowering signif- 
icance to cr = 0.0040, since the wings of the line are heavily 
blended and asymmetric. Furthermore, the core of this line is af- 
fected by non-LTE effects and was therefore removed from the 
measured line profile. 

In every line, we aimed to take full account of line blends 
and included as many data points reaching as far as possible into 
the continuum, since ZB affects the far line wings most strongly. 
However, some line wings had to be severed farther from con- 
tinuum than others, since unidentified blends were present; see, 
for instance, the far red wing of the Zeeman-sensitive line at 
6173.34 A Fe I. This line wing showed strong contributions from 
at least one blend that could not be identified in line databases. 
To circumvent this problem, iRiiedi et all (1 19971) included fake 
iron blends until the line shape was matched. We chose not to 
introduce additional fit parameters, and instead decided to sever 
the far red line wing. 



3. Analysis 

We use spectral line inversion as the tool to explore detectabil- 
ity of ZB in this paper. Spectral line inversion denotes a com- 
bination of two things: forward radiative transfer calculations 
and an iterative non-linear least squares fitting process. In other 
words, spectral line profiles are calculated using a given atmo- 
spheric model and model parameters are iteratively adjusted un- 
til the highest possible agreement between calculated and ob- 
served spectral lines is reached. The resulting best-fit model pa- 
rameters thus characterize the adopted atmospheric model and, if 
the model is assumed to be correct, the actual stellar atmosphere. 
Hence, we investigate whether introducing a magnetic field sig- 
nificantly improves the agreement between observed and calcu- 
lated line profiles. We explain the general approach of spectral 
line inversion in the following subsection. 

We further investigate model-dependence and the influence 
of model assumptions. To this end, two different kinds of models 
are used. These fitting strategies are documented in Sect. 13. 21 

3.1. Spectral line inversions 

We employ the spectral line inversion package SPINOR de- 
veloped and maintained at ETH Zurich and MPS Katlenburg- 
Lindau for this analysis. In short, SPINOR computes synthetic 
intensity profiles (SIPs), L(A), by solving the radiative transfer 
equation on a grid of optical depth points r, accounting also for 
turbulence, rotation, magnetic fields, instrumental effects, etc. It 
then determines the model parameters used to compute L(A) that 
best fit the observed data using a non-linear least squares fitting 
algorithm. We now briefly dis cuss how the S I Ps used for inver- 
sion are calculated and refer to Frutige r^t alj d2000l [2005. from 
hereon FruOO and Fru05, respectively) for further information 
regarding SPINOR. 

First, annular SIPs, I (A,0), are computed as a sum of 
weighted Voigt profiles (cf. Grav 2008, from hereon Gr08) com- 
bining thermal and pressure broadening, microturbulence (mass 
motions on a scale that is much smaller than the optical depth 
lead to a Doppler broadening due to a Gaussian velocity dis- 
persion), v m ; c , and the full Zeeman pattern as described be- 
low. "Annular" here denotes the ^-dependence, since synthetic 
profiles are calculated for each of 15 concentric annuli that 
taken together approximate a flux profile resulting from a much 
larger number of intensity profiles at different viewing angles 
(Rue di et al.lll997l Fru05). By calculating annular SIPs at dif- 
ferent angular distances from disk center, we can include limb 
darkening due to its dependence on fi = cos 6. 

To calculate t he radiation tr ansfer, we start with tabulated 
LTE atmospheres dKuruczl 19921) . These atmospheres list the pa- 
rameters temperature, pressure and electron pressure for a num- 
ber of t values and are calculated for given effective temper- 
atures, T e ff, and continuum opacity at 5000 A. In general, the 
temperature values are allowed to vary in order to obtain an op- 
timal fit to the spectral lines. Then, the whole temperature strati- 
fication T(t) is shifted up or down. At each step in the inversion 
the atmosphere is always consistently maintained in hydrostatic 
equilibrium. Hence, from a prescribed temperature stratification, 
the gas and electronic pressures, the gas density and continuum 
opacity are computed assuming hydrostatic equilibrium. Note 
that temperature values returned by SPINOR correspond to unit 
optical depth calculated for continuum opacity at 5000 A. 

The Zeeman shifts and amplitudes of the polarized c ompo- 
nents are calculated as explained in, e.g. 1 Wi ttmannl (1 1 9 74l) . A de- 
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tailed description of the implementation of the Zeeman effect in 
SPINOR can be found in Frutiger (2000) and references therein. 

In a next step, the individual annular I(A, 9) are convolved 
with radial-tangential macroturbulence profiles, ®{A, 9), that are 
intended to account for large-s cale motion d ue to convection. 
These profiles are described in iGravl d 1975b and Gr08. Thus, 
we obtain new annular SIPs, I' (A, 9) = I(A, 9) * ®(A, 9), where 
"*" denotes convolution carried out as multiplication in Fourier 
space. In accordance with Fru05, we assume equal surface frac- 
tions covered by radial and tangential flows. 

The radial velocities of individual surface elements vary 
across the stellar disk as a function of distance to the rotation 
axis, leading to rotational broadening in the integrated spectrum 
of a point source. This dependence is taken into account by intro- 
ducing 15 new angles, <p, so that the radial velocity is constant 
at any given (f>. Thus, the I' (A, 9) are turned into I' (A, 9, <f>) and 
each one is shifted for radial velocity. The disk-integrated line 
profiles are obtained by weighting the individual I' (A, 9, <p) for 
surface coverage, and summing up. 

We obtain the final SIP, L{A), by computing L(A) = V(A) * 
/i ns t, that is by convolution of I'{A) with an instrumental broad- 
ening function, /; nst . 7j nst is a Gaussian function with a FWHM 
defined as a velocity, Vi nst , in km s _1 that corresponds to the re- 
solving power of CES. Hence, we have R = 220000 = A/AA = 
c/Vi ns t. Therefore, for CES, v; nst = 1.36kms _1 . 

Summing up, our model parameters are: temperature T\ (and 
T% for two-component inversions, cf. Sect. 13.2.2b . in Kelvin at 
unit optical depth; surface gravity log g in cgs units; projected ro- 
tational (v sin /), microturbulent (v m j c ) and macroturbulent (v mac ) 
velocities, each in km s . logg was kept fixed at literature val- 
ues in most inversions, since narrow Fe I lines do not provide 
strong constraints on surface gravity. Inversions of line set B4, 
however, allowed log g to vary, since the pressure-sensitive Ca I 
line did provide this constraint. Finally, lines were allowed to 
be shifted in velocity space by < 500 m s -1 (~ 0.01 A at 6000A) 
to correct for slight differences in literature and data line center 
wavelengths, as well as possible shortcomings in the dispersion 
relation. 

Our treatment is state-of-the-art for a line inversion tech- 
nique and motivated by computational efficiency without depart- 
ing from physical realism. However, it has to be kept in mind that 
the line profiles described here are approximations that rely on a 
particular implementation (e.g. of the radial-tangential macrotur- 
bulence, or the height-independence of other parameters). This 
will be a limiting factor for the achievable fit precision and we 
can expect to find some residual signatures of this treatment in 
our analysis. Thus, an important part of interpreting our results 
will be linked to distinguishing these effects from a ZB signa- 
ture. 



To investigate detectability, we chose two different strate- 
gies. First, models with a single homogeneous atmosphere were 
considered, so called one-component (OC) models, where only 
the product Bf was determined. Second, two-component (TC) 
models were used, where B and / are naturally segregated, since 
/ defines the relative surface fractions covered by the respective 
components, cf. Eq.[3] A further refinement of this model takes 
into account the generally different temperature stratifications 
of magnetic and non-magnetic regions. This temperature differ- 
ence is well-known from the Sun (see, e.g., Solanki 1993 for an 
overview). This dichotomy in fitting strategies is addressed in 
the following and results are presented in separate subsections. 

For both strategies adopted, we allow all fit parameters aside 
from B f (OC) or B and / (TC) to vary freely in the fit (see last 
paragraph of Sect. 13. U . However, as explained above, \ogg was 
frequently used at a fixed literature value. Thus, we calculate^- 2 , 
cf. Eq.|2] for a range of values of Bf that we keep fixed in each 
inversion. In this way, we aim to avoid degeneracies that prevent 
the non-linear least-squares algorithm from finding the global 
minimum. This strategy implies that we do not aim to determine 
the other model parameters accurately. Our goal is to investigate 
whether introduction of magnetic flux (and how much) in the 
atmospheric model yields significantly better results. 

To quantify our results, we determine the overall best-fit so- 
lution, i.e. where x 2 reaches its global minimum , as well as for- 
mal errors on Bf following Pres s et all dl999b . The formal la- 
range is adopted as our standard error. To claim detection of 
magnetic fields, we require the non-magnetic best-fit solution 
to lie outside the 3<x confidence level (CL). Non-detections are 
characterized by 3cr upper limits. 

The CLs on x 2 are of a formal nature, since they depend on 
y 2 • , cf. Sect. 12.21 Errors deduced from these CLs would there- 
fore represent the true errors only in the case of a "perfect" fit. 
In addition, the choice of a specific (Kurucz) input atmosphere 
introduces a bias that affects x 1 an d thereby influences the for- 
mal error determination on Bf. However, spectral lines differ 
in their respective sensitivities to magnetic fields (usually ex- 
pressed by their effective Lande factor, cf. TableQ}. We employ 
both Zeeman sensitive and insensitive lines in simultaneous in- 
versions. Each line's different reaction to a certain postulated 
magnetic flux directly impacts the obtained best-fit^ 2 for this Bf 
value. Therefore, provided the inversions are not dominated by 
systematic effects (see, however, Sect. 14. 2. Ti l, the CLs do allow 
us to evaluate whether evidence for a magnetic field is found, 
despite the above statements about their mostly formal nature. 
The CLs further provide an order-of-magnitude estimate on the 
uncertainty involved in the measurement. 



3.2. Fitting strategies and atmospheric models 

Given the highly complex field topologies expected for sun-like 
stars, most past techniques to measure ZB assumed that coverage 
by magnetic regions was homogeneous and atmospheres identi- 
cal in all but their magnetic properties. Usually, both magnetic 
field strength, B, and surface area covered by magnetic regions 
(filling factor), /, were then determined independently, while the 
effects of multi-modal surface temperature distributions were ne- 
glected. These models computed the total radiative flux from the 
star as the sum of flux from magnetic and non-magnetic (quiet) 
surface fractions with identical temperature stratification: 

F to t — (1 - f)F qv iet + fFmagn.- (3) 



3.2.1. One-component inversions 

In the OC approach, we scan the dependence ofx 2 on the single 
parameter Bf. This is done by performing inversions for a set of 
fixed Bf values. Best-fit solutions for^ 2 and model parameters 
are obtained from each inversion. 

In OC inversions we formally set / = 1 and investigate the 
line shaping effect due to the product Bf. Therefore, we as- 
sume a mean magnetic field covering the entire surface, aiming 
to avoid the degeneracy between B and /. The spacing of the Bf 
grid is usually 50 G, and is refined to 25 G close to the overall 
best-fit x\( m for better sampling of^- 2 in this region. 
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3.2.2. Two-component inversions 

Following the solar paradigm, magnetic fields in cool stars are 
expected to be concentrated in distinct elements that are usu- 
ally d escribed as flux tubes (iStenfloll 1 973b ISpruitl[l97a : ISolankil 
Il993l) . Our two-component (TC) inversions aim to take into ac- 
count this physical fact by composing the atmospheric model of 
a magnetic and a field-free component. Thereby, we can investi- 
gate the impact of the / = 1 assumption made for OC models as 
well as the effect of a bimodal surface temperature distribution. 

Hence, we perform three cases of TC inversions: those where 
both components have identical temperature, as well as those 
where magnetic regions are either cool, or warm with respect to 
the non-magnetic component. The equal temperature case was 
frequently used in the literature and allows for a more direct 
comparison of our results than OC inversions. Our incremen- 
tal approach of increasing model complexity aims to distinguish 
between geometrical and temperature-related effects. 

We compute the best-fit x 2 for each point of the Z?-/-plane 
defined by a grid of fixed B and / values. As a result, this ap- 
proach is computationally much more expensive. 

In the TC inversions, the model contains surface filling fac- 
tor, /, as an additional fixed parameter, since we scan the B-f- 
plane. In addition to /, we introduce the additional fit parameter 
Ti for temperature at unit optical depth for the second compo- 
nent. For the equal temperature case, T2 and T\ are fitted, but 
kept equal. Both parameters can vary freely and independently 
in the cool and warm cases. To quantify our results, we construct 
formal CLs on the two-dimensional^ 2 -maps in analogy to the 
OC inversions. 

In this paper, we perform TC inversions for the Sun and 
59 Vir. For 59 Vir, for which a magnetic field is detected with 
a high CL, we investigate and compare the three different TC 
cases discussed above. 



4. Results 

We present our results in two kinds of figures, namely line fit 
plots and x 2 -plots or -maps. Figure[2]is an example of line fit 
plots that help to assess the qualitative difference between mag- 
netic and non-magnetic best-fit models. The upper panels con- 
tain the measured line profiles as ticks of length cr, overplotted 
by best-fit solutions. Solid blue represents the overall best-fit so- 
lution, and dash-dotted red lines are used to draw either the non- 
magnetic, i.e. OG (if evidence for magnetic fields is found), or 
formal 3<x upper limit on Bf. Residuals are drawn in histogram 
style in the lower panels using the identical color and line style 
theme and visualize the difference between measured and calcu- 
lated line profiles, scaled by factor 100. The purple error bar at 
the right end in each of the residual panels indicates cr, . 

The residuals contain smooth solid green lines that indicate 
the difference in line shape between the overall best-fit and the 
other best-fit model plotted, i.e. the change in line shape due 
to introduction of magnetic flux. The other fit parameters such 
as vsinj, etc., adapt to each fixed Bf value, since they are al- 
lowed to vary freely. Hence, the green line does not represent 
the Zeeman pattern alone. It illustrates the combined effect due 
to ZB and different best-fit parameters for the two cases plotted. 

Our inversions converged after relatively few iterations and 
although a given choice of starting values for the fit parameters 
(especially temperature) introduces a bias (if the x 2 surface has 
multiple minima), no dependence on these was found for the 
best-fit parameters. 



Table[3] lists our results for Bf and overall best-fit param- 
eters, i.e. those at v 2 . , for both OC and TC inversions. The 

' /Minn' 

overall best-fit values for B and / from TC inversions are found 
in the text in the respective subsections. Col. 1, 2 and 3 con- 
tain the star's name, HD number and spectral type. Spanning 
these columns, we indicate if and for what kind of TC model re- 
sults were obtained. The line sets used, as introduced in Table[T] 
are listed in Col. 4. Col. 5, 6 and 7 list the overall best-fit y 2 . , 

71 min 

number of degrees of freedom jVdof, and continuum signal-to- 
noise (S/N) measured in the data. Cols. 8, 9 and 10 provide lit- 
erature values for effective temperature, T e ff, as well as overall 
best-fit temperature at unit optical depth computed for contin- 
uum opacity at 5000 A, T\, and, in the case of TC inversions, I2, 
in Kelvin, cf. Sect. 13.21 Col. 1 1 contains literature values for log- 
arithmic gravitational acceleration in cgs units, unless marked by 
an asterisk (*). In this case, \ogg was used as a free fit parame- 
ter, constrained by the pressure-sensitive Ca I line, and the over- 
all best-fit value is presented. Col. 12, 13 and 14 list the overall 
best-fit parameters obtained for macro turbulence v mac , micro tur- 
bulence v m ; c , and projected rotational velocity, v sin i, each in km 
s . Col. 15 contains literature values for Bf, where available. 
Col. 16 lists our results for average magnetic field with formal 
lcr errors obtained from interpolation of the x 2 -plots or -maps. 
The last Col. 17 quotes the 3cr range of Bf rounded to 50 G. 

Additional figures for OC results obtained can be found in 
the online appendix. These include behavior of best-fit parame- 
ters with Bf as well as line profile plots for all OC results listed 
in Table[3]that are not included in the main article body. 

4.1. One-component inversions 

OC inversions were carried out for 4 stars in different line sets, 
see TableQ] Figure[3]shows the x 2 -plots determined for (from top 
left to bottom right) the Sun (A3), HD 68456 (A3), 61 Vir (B3), 
HD 68456 (B4), 59 Vir (B3) and 59 Vir (B4), where brackets in- 
dicate line sets used. The dashed red horizontal lines indicate the 
formal 3cr CL defined by^ n + 9, cf. Sect. l3T2~Tl 

4.1.1. The Sun 

The spectrum of sunlight reflected from Ganymede's sur- 
face is the highest-quality spectrum used in this work with 
S/N ~ 700. We a dopted sol ar abundances as determined by 
iGrevesse & Sauval d 1998I) and lKuruc3 (1 1 9921) input atmosphere 
asun. datl3. Surface gravity was fixed at \ogg = 4.44. 

The upper left panel of Fig. [3] clearly shows that the non- 
magnetic model yields our overall best-fit result for the Sun. 
We determine Bf = +90 G and a formal 3<x upper limit: Bf < 
200 G. x 2 increases smoothly and steeply above 200 G, indicat- 
ing a robust non-detection. The overall best-fit parameters are 
listed in Tableland are consistent with the literature values pub- 
lished by V&F05: vsinr = 1.63 km s -1 , v mac = 3.98 km s -1 , 
^mic = 0.85 km s _1 . 

Figure|2] shows the measured line profiles (black tick marks) 
of line set A3 overplotted by non-magnetic (blue) and magnetic 
(red, Bf - 225 G) best-fit calculated profiles. The agreement 
between either model and the measured line profiles is so close 
that amplified residuals have to be considered to view the dif- 
ference. The most Zeeman-sensitive line, CrI with g e g = 2.00 
(center), exhibits the strongest reaction to magnetic flux. This is 
indicated by the smooth green line. However, the amplitude of 
this reaction is of the same order of magnitude as <x,. 
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-0.14 +0.14 -0.14 +0.14 -0.14 +0.14 




-0.14 +0.14 -0.14 +0.14 -0.14 +0.14 
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Fig. 2: Sun, line set A3, one-component model: Data and best-fits with Bf — 0G (solid blue) and Bf = 225 G (dash-dotted red) 
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61Vir 115617 G6V B3 198.1 175 400 5571 6284 4.47 3.09 < 0.10 0.46 600 c +8 ° 0-150 



Truiillo Bueno et ID i2004h : b iLinskv et alJ dl994h ; c iGravl dl984h : * log g used as free fit parameter 



Table 3: Objects investigated sorted by increasing HD number. Values for T e g and logg from literature, unless otherwise indicated. 
Note that T\ and 7^ are SPINOR best-fit temperatures corresponding to unit optical depth calculated for continuum opacity at 5000 A. 
Overall best-fit results as well as 1 it and 3cr ranges determined for B f, see text. 3cr Bf range rounded to 50 G. 



A clear systematic signature is seen in the residuals, notably 
in the 5778.46 A Fel line. This W-shaped signature is not re- 
moved by ZB in this OC inversion and is roughly twice as large 
as the solid green line that indicates the difference between the 
non-magnetic and magnetic models. This indicates a systematic 
shortcoming in our ability to reproduce line shape and may indi- 
cate that, for instance, our approximate treatment of convection 
(as a macroturbulence) may be inadequate at this level of preci- 
sion. 

From this analysis, which is consistent with the upper limit 
of 140 G set bv lStenflo & Lindegrenl(ll977l) . also using broaden- 
ing of int ensity profiles, we can neith er confirm nor rule out the 
result bv lTruiillo Bueno et all d2004l) of Bf ~ BOG, since this 
value is consistent with our formal 3<x CL. The same is true for 



other, generally smaller, literature values of average ( turbulent) 
magnetic fields determined for the quiet Sun dSolankil2009l) . We 
conclude that magnetism averaged over the solar surface is too 
weak to be detected using this line set and method. Note that the 
above literature values refer to the quiet Sun, whereas the Sun 
was rather active at the time the analyzed data were recorded 
(13 October 2000). However, adding the active-region magnetic 
flux would raise the B f value by a couple of 10 G only and thus 
would not affect our conclusions. 



4.1.2. 61 Vir-HD 115617 

The spectrum of the G6V field dwarf 6 1 Vir in data set B was 
an ideal candidate for a robust magnetic field measurement due 
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to its slow rotation (R&S03: vsin/ < 2.2km s '), close-to solar 
metallicity ([Fe/H] = 0.05), r eff = 5571 K (latt er two in V&F05) 
and lo w X-ray activity (NEXXUS2 database dSchmitt & L iefke 
2004): log Lx — 26.65). We used an input atmosphere calculated 
for r eff = 5500 K. logg = 4.47 (V&F05) was used as a fixed 
parameter. 

The x 1 -plot in the center left panel of Fig.[3]confirms the ro- 
bust non-detection for 61 Vir. This is in agreement with expecta- 
tions moti vated by findings tha t link X-ray luminosity with mag- 
netic flux (Pevtsov e t al. 2003). The overall best-fit is reached for 
Bf = 0G and formal lcr and 3cr upper limits are established at 
Bf = 80 and 150 G, respectively. See Table[3]for details on other 
best-fit parameters. The lower 3cr upper limit obtained for 61 Vir 
than for the Sun can be explained by the higher Zeeman sensi- 
tivity of the line set B3 (61 Vir) relative to A3 (Sun). 

We found one previo us me a surem ent of magnetic flux for 
61 Vir in the literature. iGravl (11984 determined B a/Ao = 
600 G ±10%, where Ao is related to / by average viewing an- 
gle. This result is inconsistent with our null detection. We spec- 
ulate that the difference between the two measurements is due to 
our higher data quality and different, possibly more sophisticated 
analysis technique. 

4.1.3. 59 Vir - HD 11 5383 

59 Vir is a GOV field dwarf known for its activity and has 
been extensively studied by many authors. V&F05 measured 



its vsin/ = 7.4km s" 1 , T efS = 6234K, logg = 4.60 and 
[Fe/H] = 0.28 (relative to solar). logL x = 29.41 (NEXXUS2), 
approximately 2.8 dex higher than that of 61 Vir. Given its spec- 
tral type and observe d rotation period P I0t = 3.3 d measured by 
Dona hue et al.l (Il996l) . 59 Vir would be c lassified as a young star . 
However, using multiple age-indicators, Lachau me et al.l |l999) 
favor an age of 3.8 Gyr. Thus, 59 Vir is an interesting object to 
study regarding the age-activity relationship. 

Inversions for 59 Vir were performed using line sets B3 and 
B4, cf. TableQ] In the latter case, surface gravity was well- 
constrained, despite the less-than-ideal fit to the Cal line, and 
\ogg = 4.49 was obtained, consistent with the literature value. 
The other best-fit parameters are also cons istent with the litera- 
ture dTakeda et al.ll2005HSousa et al.ll2008l) . 

We determine Bf = 500+$ G for 59 Vir from line set B3 
and Bf = 525+^ G from B4, cf. the bottom left and right pan- 
els in Fig. [3] The 3<x CLs determined from both data sets are 
roughly the same with approximately 300 < Bf < 650 G. x 2 
varies smoothly over the entire Bf range, and a steep incline is 
seen above Bf ~ 800 G. Furthermore, the best-fit solution for 
Bf = G lies nearly 5cr from^ in . 

The 6173 A line was assigned a slightly higher statistical 
weight, since it is the most Zeeman-sensitive line in the set with 
geff = 2.5. The blue wing of this line contains a blend due to 
Eu II that was not fully reproduced, despite our effort to include 
all available atomic line information in the inversions. Due to the 
presence of unidentified blends, the red wing had to be severed 
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Fig. 4: 59 Vir, line set B3, one-component model: Data and best-fit with Bf = 500 G (solid blue) and Bf — G (dash-dotted red) 



before reaching continuum. Line blends as well as the systematic 
signature seen also in the results for the Sun are clearly visible. 

The residuals in Fig. [4] indicate improvement in fit quality 
due to introduction of magnetic flux for all lines. However, the 
largest differences between the magnetic {Bf = 500 G) and non- 
magnetic best-fit results are of the same magnitude as the S/N 
of the data. This result thus illustrates how difficult it is to detect 
ZB in optical Stokes /, even for relatively high values of Bf. The 
presence of systematic residual signatures similar to those found 
for the Sun and 61 Vir underline the need for rigorous and ac- 
curate implementation of radiative transfer and line broadening 
functions. 

4.1.4. HD 68456 

The F6 dwarf HD 68456 was of particular interest in our analy- 
sis, since data were available in both data sets and thus allowed 
our results to be tested for consistency. Since the dates of ob- 
servation were nearly 4 years apart, we did not simultaneously 
invert lines from both sets. After all, changes in atmospheric pa- 
rameters, for instance due to activity cycles, within this time 
span cannot be e xcluded. Fundamenta l param eters of this ob- 
ject are given by lHauck & Mermilliodl dl998l) (r eff = 6396 K, 
logg = 4.14, [Fe/H] = -0.29) and R&S03 (vsiru = 9.8km 
s -1 ). NEXXUS2 indicates a high level of activity of this star 
(log Lx = 29.05). Furthermore, HD 68456 is a particularly inter- 
esting object to study in this context, since we are not aware of 
any previous ZB measurement in a late F-type stafl 

Relatively fast rotation and high effective temperature com- 
bine to produce very shallow absorption lines in this star's spec- 
tra. This fact complicated analysis of line set A3 where lines are 
generally shallower than in data set B. Furthermore, data set A 
for this star was of lower quality than the other spectra analyzed 
in this work, containing individual data points that deviate sig- 
nificantly from the general line shape, while continuum S/N was 
approximately 400 in both data sets. The deviations seen are not 
consistent with the signature expected from a starspot. 

1 but note the work of lDonati et alj d2008f) on th e changing magnetic 
field topology on r Boo (F7 IV-V dGrav et alj|200lh ) 



The top and center right panels in Fig.[3]show our results ob- 
tained for line sets A3 and B4, respectively. For line set A3, the 
overall best-fit Bf = 1100G, but no detection can be claimed. 
This is underlined by an F-test probability Pf = 0.21 for this 
result indicating spurious fit improvement. Using line set B4, 
however, we do find evidence for a magnetic field and obtain 
Bf ~ 1 kG, with the formal 3<x CL being 600 < Bf < 1200 G. 

4.2. Two-component inversions 

The concept of two-component (TC) inversions was introduced 
in Sect. 13.2.21 We carried out these more complex inversions for 
the Sun and 59 Vir. For the Sun, for which no one-component 
detection of a magnetic field was made, we did not expect re- 
liable 2-component results. These inversions were carried out 
purely as a test. In the illustrative case discussed in the first 
subsection below we made the starting assumption for mag- 
netic flux to be concentrated in cool regions, such as spots. 
The results are quite instructive regarding systematic effects. For 
59 Vir, where the probability of attaining additional information 
from a two-component inversion is higher, we investigated the 
three cases mentioned above: equal temperatures for both com- 
ponents; warm magnetic component (warm case); cool magnetic 
component (cool case). 

4.2.1. The Sun 

Figure|5]shows the^ 2 -map calculated from TC inversions of so- 
lar line set A3. The input atmospheres used were asun.datl3 
(cf. Sect. l4. 1 .TT > for the warm and a T e ff = 4750 K Kurucz at- 
mosphere for the cool component. We obtain = 170, a 
value considerably lower than that obtained for the OC model, 
for B = 4950 G and / = 0.05. The Bf grid was calculated for 
field strengths of up to 5000 G, with filling factors up to 0.95 
for B < 2000 G and / < 0.5 for B > 2000 G. The CLs are 
235 < Bf < 250G(lcr) and 150 < Bf < 1400G(3cr), and 
thus would suggest evidence for very strong surface fields in the 
Sun. However, there are a number of problems with these re- 
sults. For instance, magnetic fields in sunspots do not usually 
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exceed ~ 3kG while th e average field strength o f sunspots is 
roughly 1000 - 1300 G (ISolanki & Schmidtll 19931) and the field 
strength averaged over sunspot umbrae usually does not exceed 
2500 G. In contrast, our 3cr confidence level suggests that B is 
well above 4kG with / between 1— ~ 30% in disagreement 
with the observed sunspot coverage which is (always) below 
1%. Furthermore, a multitude of measurements places the av - 
erage Bf of the (quiet) Sun considerably lower (Solanki 2009). 
Last, the 3cr range for Bf is vastly larger than that found for 
the OC result. Therefore, we conclude that our fit was not well- 
constrained and that in this case the more complex TC model 
yielded less significant results than the simpler OC model. 

Figure[6]helps to understand what happened here. The great- 
est decrease in^ 2 was achieved by decreasing the systematic sig- 
nature in the 5778 A Fe I line. This apparent improvement in fit 
quality has stronger influence on^- 2 than the deteriorating agree- 
ment between the measured and calculated profiles of the CrI 
line at 5783 A, since the former line contributes more degrees of 
freedom. Therefore, the magnetic field (or rather the lower tem- 
perature associated with the magnetic component) improves the 
fit of the Zeeman-insensitive (g e g = 1 .2) line, while the Zeeman- 
sensitive (g e ff = 2.0) line is less well matched. The worse fit in 
the blue Cr line wing is compensated in part by a small unidenti- 
fied blend in the red wing that is falsely removed by the magnetic 
model. Thus, our conclusion is that this TC inversion is driven 
by systematic difficulties and yields a result that is less repre- 
sentative of the physical processes, albeit superior in terms of 
pure^ 2 . It therefore serves as an example of how line blends and 
other systematic effects on line shape affect our results. 

4.2.2. 59Vir-HD 115383 

The results obtained for 59 Vir using OC models indicated for- 
mal 3<t detections of significant magnetic flux. To assess, if this 
detection is actually robust and significant, we performed TC in- 
versions for this strongest candidate. The following results were 
obtained from our physically motivated TC inversions that dis- 
tinguish between magnetic and non-magnetic components. 

We investigated three different cases of TC inversions: equal 
temperatures, cool magnetic, and warm magnetic components. 
Thus, in the latter two cases, we investigate the effect of having 
different temperatures for the two components. 

For the equal temperature case, we used the T e g = 6250 K 
Kurucz atmosphere for both components and coupled the two 
parameters T\ and Ti, forcing them to be equal. In the cool case, 
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Fig. 5: ^ 2 -map for the Sun, line set A3, TC model with cool mag- 
netic regions. 



we used tabulated atmospheres calculated for = 6250 K and 
5000 K for the non-magnetic and the magnetic component, re- 
spectively. The warm case was motivated by the warmer plage 
or network regions on the Sun. The Kurucz atmospheres were 
those with T e ff = 6250 K for the magnetic and T e ff = 6000 K for 
the non-magnetic regions. The temperature parameters T\ and 
Ti could vary freely and independently in both the warm and 
cool cases. The fitting ranges for T\ and were large enough 
to allow turning the cool component into the hotter one and vice 
versa. The starting values of the fit parameters were identical in 
both cases, apart from the inverted values of T\ and T2 (the sec- 
ond temperature component was magnetic). All other starting 
values were identical to the OC inversions of 59 Vir. 

Equal temperatures: FigurePTal shows the ^ 2 -map obtained 
from our inversions of a model with equal temperatures for both 
components. As can be seen in Table|3] our overall best-fit pa- 
rameters are virtually identical with our OC results. Even the two 
B f values lie within their lcr CLs. v 2 . is reached for B = 700 G 

J min 

and / = 0.6. However, even the lcr CL spans a large range in 
/-values, and the 3<x CL is consistent with almost any / value. 
The reason for this lies in the x 2 surface, which has the well- 
known banana-shape, with a minimum running from the combi- 
nation of high / and small B to the opposite extreme of small 
/ and large B, while maintaining roughly equal magnetic flux 
(i.e. Bf). This is the signature of degeneracy between these two 
parameters. Of particular importance is the fact that in spite of 
the additional parameter /, the x\m is identical to the OC case, 
which also indicates that there is not sufficient information in the 
data to distinguish between B and /. Furthermore, the formal Bf 
error margins for this case are larger than for the OC inversions 
and the other best-fit parameters nearly identical. Therefore, we 
conclude that a distinction made between magnetic and non- 
magnetic regions is insignificant. 

iLinskv et alj d!994l hereafte r Li94) det ermined Bf = 190 G 
for 59 Vir (number published in lSaar| [l996) using an equal tem- 
perature TC model. Their Bf value is lower than, yet consistent 
with, our 3<x Bf range obtained from the TC inversions in the 
equal temperature case. It is, however, inconsistent with our 3<x 
Bf range determined from the OC inversions that favor higher 
Bf. At present it is not possible to say whether the difference to 
our results is due to differences in technique or due to the evolu- 
tion of the stellar magnetic field. 

Cool magnetic regions: As in the above example for the Sun, 
we assume the majority of the magnetic flux to originate in 
starspots, i.e. in cool regions on the stellar surface. FigurelTbl 
shows the ^- 2 -map corresponding to this case. The overall best- 
fit Y 2 . = 216 is reached for B = 250 - 300 G with f = 0.475 

xv min ■> 

for both, i.e. Bf - 120 - 140 G. x^ is lowered by 46.3 relative 
to the OC and TC equal temperature models and the formal 3cr 
confidence level is < Bf < 320 G. The formal 3<x range in 
Bf is smaller than in the equal temperature TC case. The im- 
provement in fit quality over the OC model is most clearly seen 
in the better matching central depth of the Zeeman-insensitive 
6151 A and of the EuII blend to the Zeeman-sensitive line at 
6173 A, cf. Fig. [8] Both improvements are likely due to changes 
in level populations. The presence of a cool component enables a 
better matched line depth of Fe I, whereas the warm component 
contributes additional ionized europium thereby better matching 
the Eu II blend. In fact, the other two iron lines are also better 
matched by the TC inversion. 
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Fig. 6: Sun, line set A3, comparison of overall best-fit TC (B = 4900 G, / = 0.05, blue solid) and OC (Bf = G, red dash-dotted) 
results. 





500 1000 1500 2000 500 1000 1500 2000 
B [G] B [G] 

(a) ^ 2 -map, same atmospheres for magn. and (b)^ 2 -map, cool magnetic regions 
non-magn. regions 




500 1000 1500 2000 
B[G] 

(c)^f 2 -map, warm magnetic regions 



Fig. 7: Results obtained from TC inversions of 59 Vir, line set B3. x\«& indicated by large cyan pluses. CLs indicated by solid white 
lines. 



Figure|9] clearly shows that the difference in best-fit Ti and 
T 2 within the 3<x CL lies between 1400- 1600 K. This is roughly 
in agreement with the difference in temperatures chosen for the 
input atmospheres and the temperature differences between the 
quiet solar photosphere and sunspot umbrae. 



Warm magnetic regions: For this case, we employed a model 
containing slightly warmer magnetic than non-magnetic regions, 
i.e. assumed a model representative of plage or network-regions. 
This case is particularly interesting, since a warm second com- 
ponent contributes more flux than a cool one. Thus, the effect 
due to a certain B f on line shape can be expected to be stronger. 

Fi gurelTcl illustrates our results. = 226 is slightly higher 
than in the above cool case and yields B = 450 G and / = 0.6. 
Bf ranges from 210 - 325 G in the formal lcr interval, while we 
find it consistent with - 450 G at 3<x. 



The 3cr range in B values is roughly the same as in the cool 
case, whereas the 5<x level is limited to lower B values. Opposite 
to the cool case, high-/ solutions are favored. Thus, both cases 
favor a slight majority of the surface covered by warm elements. 
Despite the small (250 K) difference in T e g of the input atmo- 
spheres, the difference in best-fit parameters T 2 and T\ within 
the 3cr CL. ranges 1300 - 1400 K (AT = T 2 -T\ in the warm 
case, since T 2 > T{), This temperature difference is very similar 
to that found in the case of cool magnetic regions (for the 3<x 
CL). 



Comparison of TC cases (equal, cool & warm): The equal 
temperature TC case yields results in agreement with OC in- 
versions, but wider CLs on Bf. Significantly lower x*. is found 
in the two different temperature cases, also for the non-magnetic 
best-fit solutions (lowest^ 2 for Bf — 0G in equal, cool & warm 
cases: 287, 217, 231). Interestingly, the warm and cool TC cases 
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Fig. 8: 59 Vir, line set B3: Comparing Bf - 120 G TC (blue solid) to Bf = 500 G OC (red dash-dotted) best-fit results. 

5. Discussion 

To assess detectability of ZB for this kind of analysis, we first 
discuss the general uncertainties and sources of error before ex- 
amining the impact of model assumptions. 

Our analysis is prone to uncertainties from various sources, 
such as degeneracies among various line broadening agents and 
uncertainties in information used as input (atmospheres, line 
data, lines used, etc.). We attempted to avoid degeneracy among 
line broadening agents by allowing all fit parameters to vary 
freely for given average magnetic fields. This approach helped 
to nudge the non-linear ^-minimization routine towards the 
global minimum. In addition, it provided a clear visualization 
(see PigEb °f the range of acceptable Bf values, thereby pro- 
viding very solid upper limits on Bf of, say, a few hundred to 
1 .5 kG depending on the star. 

More spectral lines could have provided more stringent con- 
straints on model parameters. Unfortunately, our choice of lines 
was severely limited by the small wavelength coverage of the 
CES spectra. Temperature stratification would have been even 
better constrained by using lines from multiple ionization stages. 
In order to achieve the fit quality presented in this paper, how- 
ever, we relied on lines from different excitation potentials. 
We had found that some combinations of lines yielded better 
fit agreement than others. The reasons for this empirical fact 
are likely related to blends; specifically, to blend abundances 
(often unconstrained), lack of literature on transition param- 
eters needed for SPINOR, uncertainties in oscillator strengths 
or rest-wavelengths, and unidentified line blends or even rela- 
tive displacements of blend wavelengths due to convection. For 
a literature co mparison of the fit quality presented here, see 
iFrutiger et al.l d2005l) who applied SPINOR in simultaneous in- 
versions of 14 lines for a Cen A and B, including lines from 
different ionization stages. Their Figs. 2 & 3 clearly illustrate the 
larger discrepancy between model and data when a large number 
of lines is used. 

Small systematic effects (best visible in the spectral differ- 
ence between observed and computed lines profiles) compli- 
cated our analysis. The high quality of the data revealed these 
systematic shortcomings in reproducing real line profiles. The 
data for the Sun was both the highest quality spectrum (in terms 
of S/N) and one of the least rotationally broadened. The fact that 
the TC inversions for the Sun produced such implausible results 



tend towards virtually identical atmospheric configurations with 
nearly the same overall best-fit parameters for T\, T2 (though 
T\ and T2 are switched among the cool and warm cases), v sin 2, 
v mac , and v mic , cf. Table[3] 

This suggests that the true reason for the improved fit to 
the observations is the presence of two temperature components. 
These simulate the inhomogeneity of the stellar atmosphere pro- 
duced by starspots, granulation, oscillations, etc. and possibly 
also non-LTE effects in the line formation, which are not taken 
into account in the LTE modeling we have carried out. 

A key difference between the different TC cases is the flux 
contribution from the respective magnetic component. Due to 
this difference in contrast, different Bf values may be expected. 
Our inversions yield comparable Bf for the OC and TC equal 
temperature model, while B f is lower for the cases with different 
temperatures. Among these latter two, the warm case (brighter 
magnetic component) yields a slightly larger Bf. However, both 
are consistent with no magnetic field at the 3<x CL (though the 
warm case is not consistent with Bf = 0G at 2<x). 




500 1000 1500 2000 
B[G] 



Fig. 9: Ar-map for 59 Vir, cool case. Difference in K between T\ 
and T2 indicated by colors indicated in legend on right. Contours 
drawn for AT = T\ — T2 = OK (orange), 500 K (magenta), 
1400 K (blue), 1600 K (purple), and 2000 K (black). White con- 
tours indicate 3 and 5cr CLs. The large white plus indicates the 
point of overall best fit. 
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may indicate that a limit was surpassed for an adequate repre- 
sentation of line profiles by simple broadening functions such 
as the radial-tangential macroturbulence profile. However, some 
of the difficulties in the solar TC inversion were also related to 
the lower Zeeman sensitivity of line set A3, in addition to the 
low level of solar magnetism and the data-related problems of 
line blends, etc. Further candidates for the origin of the system- 
atic signatures at this level may be related to any of a num- 
ber of sources, such as the implementation of limb-darkening 
and disk-integration, temperature stratification, neglect of stel- 
lar convection (granulation) and oscillations, differential rota- 
tion, neglect of non-LTE effects, depth-dependent magnetic field 
strength, scattered light in CES data, etc. In addition, the temper- 
ature gradients of active stars may differ significantly from those 
of inactive stars. Thus, the representation of magnetic regions 
by "quiet" atmospheres may also be too crude an approxima- 
tion, since it neglects inhibition of convection and the influence 
of magnetic fields on temperature stratification. Thus, applica- 
bility of LTE Kurucz input atmospheres could be questioned. 
Non-LTE effects might also be seen, since the systematic resid- 
ual signature appears to be enhanced in hotter stars. The analysis 
thus suffers from a multitude of highly degenerate effects orig- 
inating in approximations. Most of these flaws cannot be easily 
resolved, however, and there are valid reasons for adopting at 
least some of these approximations. 

Two kinds of fitting strategies were explored: one- 
component (OC) and two-component (TC) models, cf. Sect. 13. 21 
for details. OC models assume a homogeneous atmosphere with 
a single temperature stratification and an average magnetic field, 
Bf, covering the entire stellar surface. For the Sun and 61 Vir, 
both stars with low X-ray activity, our OC results excluded Bf 
exceeding a few hundred Gauss. OC inversions further yielded 
evidence of 500 G to 1 kG fields for two faster rotators with 
higher X-ray luminosity, namely 59 Vir and HD 68456. As could 
be expected from its larger maximum g e s value and g e s cover- 
age, line set B3 was more sensitive to ZB than A3. 

Two-component (TC, cf. Sect. 13. 2. 2b inversions were per- 
formed to distinguish between magnetic and non-magnetic sur- 
face regions. For our solar spectrum in line set A3, we present a 
case study that assumed the magnetic fields to be concentrated 
in spots, i.e. in significantly cooler regions. However, the solar 
TC model was driven by systematic effects (partly due to weak 
blends in the line wings) and yielded a result inconsistent with 
the OC result and the literature. 

For 59 Vir, we explored three different cases of TC inver- 
sions: a) equal temperatures for both components; b) a cool mag- 
netic, and c) a warm magnetic component. We find that the equal 
temperature TC case is consistent with the results from the OC 
inversions. Furthermore, the Bf determined is consistent with 
the only available literature value. However, no improvement in 
X 2 is seen over the OC model, while an additional fit parame- 
ter, /, was introduced. This is consistent with the well known 
result that the product Bf i s much better constrained than B and 
/ individually dSaarlll990al) . 

The picture changes, however, when the model takes into 
account two different temperatures in the two components, as 
was done for the cool (spots) and warm (plage and network re- 
gions) TC cases. In fact, a significant improvement of x 2 was 
seen for both of these more complex inversions. Independent of 
the input atmospheres used for the two components, these in- 
versions converged towards a very similar photospheric model 
featuring similar coverage (f ss 50%) by hot and cool features 
with AT * BOOK. The difference in Ax 2 min ~ 10 between the 
warm and cool case is probably the result of the a priori choice 



of input atmospheres. Magnetic flux did not play a deciding role 
in the inversions' ability to converge towards this common over- 
all best-fit model. Note, however, that the improvement of the 
fit comes at the expense of additional free parameters (2 new 
parameters are introduced compared to the OC case). 

For both the warm and the cool TC case, Bf was lower than 
in the OC or TC equal temperature case and consistent with no 
magnetic field at the 3<x CL (however, not at 2cr in the warm 
case). Thus, while some evidence for ZB was found in these 
cases, no detection can be claimed when allowing for two differ- 
ent temperature components, one magnetic, the other field free. 

Given the significant improvement of ^^ in , we find that 
model assumptions impact our ability to reproduce the observed 
line shape more strongly than the presence or absence of ZB in 
optical Stokes /, even for relatively high Lande factors. At the 
same time, the presence of ZB is not excluded, either. 

As a next step, even more complex models might be invoked. 
However, this does not appear to be expedient at this point, given 
the limited number of lines used for analysis, the approximative 
treatment of line profiles (e.g. the radial-tangential macroturbu- 
lence and the height-independent microturbulence), the small- 
ness of the Zeeman signature, and the fact that the TC inversions 
of 59 Vir already reached the noise limit. 

A more prudent approach would be to increase sensitivity to 
ZB. This could be done by using high-resolution infrared spec- 
tra. In this way, the A 2 dependence of Eq.Q] could be exploited 
to provide stringent constraints that are in addition to the differ- 
ential ^/-response that was at the base for this investigation. 
Furthermore, the IR contrast between cool and warm regions 
would be lower, thereby enhancing Zeeman signatures (e.g. 
iJohns-Kmlll |2007|) . In M-dwarfs, successful measurements of 
strong magnetic fields have been reported using mole cular lines, 
e.g. using F eH lines in the Wing-Ford bands dReiners & Basri 
2006, 2007). With the majority of their flux in the infrared and 
strong fields, these stars constitute prime candidates for direct 
Bf measurements using spectral line inversion as outlined in this 
work. 

6. Summary 

In this paper, we investigated detectability of ZB in optical 
Stokes / spectra of the Sun and three slowly rotating sun-like 
stars by performing spectral line inversion with SPIN0R. 

Thanks to the ability of SPINOR to perform spectral line in- 
version, different scenarios could be compared and model as- 
sumptions evaluated by scanning the x 2 plane for fixed Bf, or 
combinations of B and /, values. This approach enabled us to 
investigate the form of x 2 as a function of Bf as shown in Fig. [3] 
and therefrom determine formal errors on Bf. 

Systematic uncertainties in our analysis are in part related 
to literature used as input, such as atomic line data and ele- 
mental abundances. Due to the small magnitude of the Zeeman 
broadening even weak line blends complicated the analysis. 
Approximations, such as the representation of turbulence in 
the stellar atmosphere by radial-tangential macroturbulence and 
height-independent microturbulence, also influenced our results. 

OC inversions excluded the presence of significant magnetic 
flux {Bf < 150 G at the 3cr CL.) for the Sun and 61 Vir. Our 
results for HD 68456 may be interpreted as first evidence for a 
magnetic field in a late F-type dwarf obtained by the ZB tech- 
nique, with Bf x 1 kG and formal 3cr errors of approximately 
350 G. The OC result for 59 Vir is B f = 500^° G. 

The TC equal temperature case yielded Bf = 420 G for 
59 Vir, a higher value than that reported by Li94, but consistent 
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within the (large) formal lcr error. TC inversions with compo- 
nents of different temperatures were consistent with B f between 
G and 450 G at the 3cr CL. The case of a warm magnetic com- 
ponent was inconsistent with an absence of magnetic field at 2<x. 
Previous investigations of this kind frequently used models sim- 
ilar to our equal temperature TC case. Our investigation of the 
different TC inversion cases demonstrates the strong influence of 
model assumptions. In fact, the freedom to freely and indepen- 
dently vary both atmospheric temperatures in TC models was 
shown to influence line shape more strongly than ZB for the 
lines used in this work, cf. Sect. 14.2.21 The improvement of x 2 
due to a second component of different temperature was signif- 
icantly larger than that due to a magnetic field, irrespective of 
whether the magnetic component was defined as warm, or cool. 
The present data quality exceeds that of previous studies, while 
our sensitivity to Zeeman broadening as expressed by the g e ff 
values in our line sets is comparable to that in the literature with 
the exception of studies in the infrared. 

We conclude that measurements of ZB in optical Stokes / 
data of slowly rotating sun-like stars are subject to large un- 
certainties. These are mostly due to data related uncertainties, 
assumptions on atmospheric models, approximations in line 
broadening agents, and degeneracy between line broadening 
agents. However, the Zeeman sensitivity of every individual line 
(as expressed by its value) differs from line to line by a factor 
of up to 2.5. This provides a powerful discriminant for our anal- 
ysis, since the other relevant broadening effects depend mostly 
on temperature or wavelength and are thus essentially the same 
for all lines used. It is this differential reaction of each spectral 
line to a given magnetic field that ensured the significance of the 
measurements presented in this work, despite the small magni- 
tude of the ZB signature. 

The most promising way to improve the present approach is 
to further exploit the A 2 -dependence of ZB and use longer wave- 
length data. In the (near) infrared, Zeeman broadening clearly 
dominates over Do ppler effects and Zeeman splitting can b e 
seen directly, see e.g . Saar jM994 : I.Tohns-Krull & Valentil(ll996l) : 
IValenti & Johns-Krulll (12001). New high-resolution near infrared 
spectrographs have become available for these tasks. 
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Appendix A: Online figures 

We present additional figures illustrating our OC results in this 
appendix. 

A. 1 . Additional one-component fit results 

Figure lATI shows the OC fit results not presented in the main ar- 
ticle body whose corresponding^- 2 -plots were included in Fig. [3] 
The best-fit models drawn are mentioned in the respective cap- 
tions of the sub-figures. 

A.2. Free fit parameters 

Figure lA~2l shows the behavior of the free fit parameters with 
Bf in the OC inversions. Red error bars indicate formal singular 
value decomposition (SVD) errors obtained from SPINOR. They 
also indicate the spacing of the Bf grid scanned, thereby illus- 
trating the higher sampling close to x^ nin - Generally, fit param- 
eters behave smoothly, apart from regions where small numer- 
ical jumps are visible, such as for vsim in Fig. IA.2al The cen- 
ter panel in each sub-figure labeled v tur b contains both macro- 
and microturbulent velocities, with macroturbulence generally 
stronger than microturbulence. Degeneracy of fit parameters is 
evident from the combined reaction of the parameter set caused 
by the fixed values of Bf. 
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Fig. A.l: Additional comparison of overall best-fit results from one-component inversions with non-magnetic models and 3<x upper 
limits, respectively. 
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Fig. A. 2: Behavior of best-fit parameters with Bf using one-component models. Plot organized the same way as Fig. [3] v tur b indicates 
micro- and macroturbulent velocities, with macroturbulence usually larger than microturbulence. 



